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Topological quantum error correction codes are extremely practical, typically requiring only a 2-D 
lattice of qubits with tunable nearest neighbor interactions yet tolerating high physical error rates 
p. It is computationally expensive to simulate the performance of such codes at low p, yet this is 
a regime we wish to study as low physical error rates lead to low qubit overhead. We present a 
very general method of analytically estimating the low p performance of the most promising class 
of topological codes. Our method can handle arbitrary periodic quantum circuits implementing the 
error detection associated with this class of codes, and arbitrary Pauli error models for each type of 
quantum gate. Our analytic expressions take only seconds to obtain, versus hundreds of hours to 
perform equivalent low p simulations. 



A number of proposals exist for a 2-D array of qubits 
with tunable nearest neighbor interactions [IH1]. Any 2- 
D topological quantum error correction code [SHH] and 
some 3-D codes (TUJ [TTJ can be mapped with little or no 
overhead to such hardware. Despite the undemanding 
physical requirements, most forms of topological quan- 
tum error correction (TQEC) possess very high thresh- 
old error rates. The threshold error rate of a quantum 
error correction code is the physical error rate p t h be- 
low which it is possible to perform arbitrarily reliable 
quantum computation. The best TQEC schemes possess 
threshold error rates of order 1% [12j [13] . 

Viewing the 2-D array of qubits together with the du- 
ration of a quantum computation as a 3-D volume, we 
define a TQEC circuit to be any 3-D periodic quantum 
circuit that detects errors (using combinations of single- 
qubit measurement results) in such a way that detection 
events can be unambiguously divided into two classes, 
primal and dual, and a single physical error away from 
boundaries of the 3-D volume always leads to or 2 local 
primal and dual detection events. Since distinct detec- 
tion events are associated with distinct sets of qubit mea- 
surements, and measurements begin at specific space- 
time locations, every detection event can be systemat- 
ically associated with a unique point in space-time. We 
can use our tool Autotune [2] to efficiently determine 
the total probability of single errors leading to a pair of 
(primal or dual) detection events at any chosen pair of 
space-time locations. We visualize each such probabil- 
ity as a cylinder between the relevant space-time loca- 
tions with diameter proportional to the probability. We 
call each cylinder a stick, and call any large collection of 
sticks a nest. Fig. [I] contains an example of a (primal) 
nest generated by Autotune. This nest is associated with 
the surface code [T5HT7] . Only a basic understanding of 
the definitions of detection events and nests contained in 
this paragraph is required to read this work. 

The distance d of a TQEC instance can be deter- 
mined from its nests by counting the minimum number 
of sticks required to connect distinct boundaries. Bound- 



aries in fig.[TJare indicated by sticks apparently leading to 
nowhere, so it can be seen that this example has distance 
d = 4. By analyzing a particular nest in the limit of large 
d, we have been able to prove [TB] that the surface code 
possesses a non-zero threshold error rate when minimum 
weight perfect matching [T9lf2T| is used to correct errors. 
In this work we extend of the techniques of |18j much 
further, and describe how to automatically calculate two 
analytic expressions p^ 1 (d, p) and p^ 1 (d,p), valid for even 
d, for the asymptotic low p performance of any TQEC 
circuit as defined above with any Pauli error model for 
each type of gate. Our new method requires seconds of 
computation versus typically hundreds of hours to do an 
equivalent simulation based low p analysis of a TQEC 
circuit. 

We begin with a few definitions. A logical operator 
is any path of sticks connecting distinct boundaries. A 
fault is any selection of d/2 sticks from a logical operator. 
The probability of a logical operator p op or fault pi is the 
product of probabilities of its constituent sticks. We de- 
fine the compliment probability p c to be p op /pf, namely 
the product of the probabilities of the sticks in the logical 
operator that are not in the fault. Note that more than 
one logical operator can contain the same fault. Mini- 
mum weight perfect matching based TQEC fails when 
a fault occurs that is contained in a logical operator 
such that p c > pf. The matching algorithm will insert 
corrections corresponding to the most probable selection 
of sticks that generates the observed detection events, 
which will result in a logical operator that corrupts the 
stored quantum data. Failure will also occur 50% of the 
time when p c — pf as this case corresponds to two differ- 
ent equally probable complementary faults that lead to 
the same detection events. Given the identical detection 
events for both faults, matching will randomly succeed 
50% of the time. 

Each layer of boundary sticks in fig. [TJ corresponds to 
the execution of a single round of error detection circuitry 
- a complete period of the circuit. Our goal is to for- 
mally define and determine the total probability of all 
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FIG. 1: An example of a surface code nest. Each cylinder 
represents the total probability of a single error leading to 
detection events at the endpoints of the cylinder. Time runs 
vertically. 



faults causing failure per round of error detection. 

The first step towards our goal is to imagine we have a 
nest of infinite temporal extent and to choose a particular 
round of error detection to focus on. More specifically, 
we wish to focus on the boundary sticks associated with 
a single boundary in that round. In practice, we generate 
nests with 2d + 4 rounds, as shown in fig. [I] and choose 
round d + 2 — this ensures that everything we consider 
is well away from temporal boundaries. 

Armed with our chosen boundary sticks, we sequen- 
tially generate every d-stick logical operator starting on 
each of these boundaries sticks and sequentially generate 
every fault contained in each logical operator. Every time 
we generate a new fault, we add it to a set, including its 
value of p c relative to the current logical operator. If the 
fault is already present, p c is adjusted to the maximum of 
the old and new values. This is done to make sure that 
only the most likely fault compliments are considered. 
When complete, we will have a base set of faults. 



TABLE I: Surface code distance d, stochastic simulation de- 
termined primal and dual asymptotic coefficients (Ap d ^ 2 ), 
time required to obtain the lowest p data points used to cal- 
culate the asymptotic coefficients, analytic primal and dual 
asymptotic coefficients (Bp d ^ 2 ), time required to perform the 
analytic computation. 



We need a set of faults uniquely associated with each 
round of error detection. Our base set is too large, con- 
taining many faults that should be associated with ear- 
lier or later rounds. We therefore consider the boundary 
sticks one round higher, round d + 3, and again sequen- 
tially generate every c£-stick logical operator and fault. If 
a fault is not found in our base set, we add it to a step set, 
adjusting p c if required as described above if the fault is 
already in step. If a fault is found in our base set and the 
new p c is higher, the fault is moved to the step set with 
the new p c . The step set constructed in this manner does 
not contain faults that should be associated with rounds 
in the past. It will, however, contain faults that should 
be associated with rounds in the future. 

We continue considering successively higher rounds of 
boundary sticks, sequentially generating all logical oper- 
ators and faults until we fail to generate a single fault 
matching a fault in step but with higher p c . When this 
process is complete, our step set will contain only faults 
best associated (highest p c ) with logical operators mak- 
ing use of one of the g?+3 round boundary sticks. We can 
now sum the fault probabilities pt , using a scale factor of 
if p { > pc, 0.5 if pf = pc, or 1 if p { < p c . 

Having described our method, we turn our attention to 
its performance. Fig. [2] show simulation data for a stan- 
dard surface code making use of initialization to |0) and 
|+) and measurement in the Z and X bases each with 
probability p of failure, CNOT and identity gates with 
depolarizing noise of probability p, and all gates of equal 
duration. These seemingly innocuous curves took over 
1500 CPU hours to generate, with the vast majority of 
the time spent generating the lowest p data points. De- 
spite the substantial computational effort, it can be seen 
that the d = 8 curves in particular have not yet well con- 
verged to their expected low p quartic asymptotic form 
(quartic due to any fault leading to failure involving a 
minimum of 4 sticks). The coefficients A of each simu- 
lation determined Ap d / 2 can be found in Table [TJ along 
with the time required to obtain the lowest p data point. 
Note that our simulations include the full depolarizing 
channel and hence all of our simulations return both the 
primal and dual performance. 

Table [T] also contains analytic coefficients B obtained 
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FIG. 2: Probability of (top) logical Z (primal) and (bottom) 
logical X (dual) error per round of error correction for vari- 
ous surface code distances d and physical error rates p. The 
asymptotic curves (dashed lines) are quadratic, cubic and 
quartic. 

using our new method. For distances 4 and 6, where de- 
cent asymptotics could be obtained using stochastic sim- 
ulations, the agreement is good. Discrepancies may be 
due to the assumption in the analytic calculation of in- 
dependent stick probabilities, which are not independent 
in simulations. Given a Pauli error model is a significant 
approximation of the physics of real devices, the maxi- 
mum 12% discrepancy for distance 4 dual asymptotics is 
of no practical concern. Of great practical benefit is the 
ability to analytically estimate the asymptotics of dis- 
tances out of reach of simulation. It can be seen from 
fig. [2] that the distance 8 asymptotic curve is expected to 
be lower, an expectation that is confirmed by the analyt- 
ics in table [T] We can now claim a solid understanding of 
the error processes leading to failure in TQEC schemes 
of the type considered in this work. 

From [18], we expect exponential suppression of log- 
ical error with increasing d at fixed p below threshold. 
This expectation is supported by the ratio of succes- 
sive pairs of primal B coefficients being approximately 37 
and that of the dual B coefficients being approximately 
32. This allows us to write p p L r (d,p) = 0.21(37p) rf / 2 and 
pf l {d,p) = 0.23(32p) d / 2 , valid for p < 1CT 4 . Note that 
only the distance 4 and 6 analytic data are required to 
construct these expressions, requiring 2.5 seconds of com- 



putation versus 750 hours to obtain the same data by 
simulation, a factor of a million improvement. 

In summary, we have described a highly computation- 
ally efficient method of obtaining accurate analytic ex- 
pressions describing the low p performance of topologi- 
cal quantum error correction, the form of quantum error 
correction that leads to the lowest space and time over- 
head for large quantum computations and the parameter 
region that was previously computationally almost inac- 
cessible. Our method can cope with the details of real 
hardware, including geometric constraints, gate duration 
asymmetries, gate error rate asymmetries and arbitrary 
Pauli error models for each type of gate as our method 
builds on our tool Autotune. We expect our method 
to facilitate the real-time interactive investigation of the 
dependence of the error correction performance on the 
details of the hardware, something impossible via sim- 
ulation, greatly assisting the search for less experimen- 
tally challenging pathways towards large-scale commer- 
cially viable quantum computation. 
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